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PREFACE 


This  Memorandum  is  one  in  a  continuing  series  of  RAND 
publications  dealing  with  theoretical  and  computational 
questions  which  have  arisen  in  connection  with  the  RAND 
program  of  research  in  biology  and  physiology. 

The  Memorandum  contributes  to  our  ability  to  apply 
electronic  computer  techniques  in  the  analysis  of  complex 
chemical  systems,  which,  until  recently,  were  far  too 
complex  for  any  quantitative  analysis--e .g. ,  models  of 
complex  physiological  systems.  But,  as  the  capability  to 
use  computers  in  such  applications  becomes  better  under¬ 
stood,  certain  hitherto  unresolved  questions  of  mathematics 
and  chemistry  become  apparent.  This  Memorandum  considers 
the  resolution  of  several  of  these  problems,  and  should  be 
of  interest  to  both  mathematicians  and  chemiists. 

The  mathematical  aspects  should  also  be  of  interest 
in  other  fields  where  computational  analyses  of  complex 
chemical  systems  are  under  consideration;  e.g. ,  in  studies' 
of  rocket  propulsion  systems,  planetary  atmospheres,  re¬ 
entry  problems,  etc. 

^For  example,  see  Refs.  4,  6-7,  9-12,  and  14. 
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SUMMARY 


Recent  progres.s  in  computational  techniques  for  the 
analysis  of  complex  chemical  systems  has  renewed  interest 
in  the  relationship  between  the  mass  action  laws,  which 
have  been  traditionally  employed  for  the  computacion  of 
chemical  equilibrium  compositions,  and  the  Gibbs  free 
energy  function,  which  until  recently  has  generally  been 
used  only  for  theoretical  purposes.  Several  attendant 
mathematical  questions,  hitherto  unresolved,  have  now  be¬ 
come  important.  These  concern  the  existence  and  unique¬ 
ness  of  a  solution  to  the  equations  arising  from  the  mass 
action  laws;  the  existence  and  uniqueness  of  a  composition 
that  minimizes  the  free  energy;  and  the  precise  relation¬ 
ship  between  the  two. 

These  questions  become  imp'"  . cant  because  there  is  a 
greater  probability  of  complex  pathologies  being  exliibited 
when  extremely  large  chemical  systems  are  analyzed,  and 
because  the  use  of  automatic  digital  computation  demands 
a  precise  knowledge  of  these  pathologies. 

In  Part  i,  the  problems  are  formulated,  both  from  the 
mathematical  and  chemical  viewpoints.  In  Part  II,  mathe¬ 
matical  tools  unavailable  in  Gibbs'  time  le  employed  in 
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an  attempt  to  answer  the  questions  posed.  Some  applica¬ 
tions  of  the  results  are  given  in  the  concluding  section 
for  illustrative  purposes. 
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INTRODUCTION 


The  problems  discussed  in  this  Memorandum  are  moti¬ 
vated  by  problems  arising  in  the  computation  of  chemical 
equilibrium.  Classically,  the  composition  of  a  single-  or 
multi-phase  chemical  system  in  chemical  equilibrium  has 
been  computationally  determined  by  solving  certain  systems 
of  simultaneous  equations.  These  equations  included  linear 
equations  (mass  balance  laws)  as  well  as  nonlinear  equa¬ 
tions  (mass  action  laws).  Gibbs  [1]  and  others  formulated 
the  "free  energy  function"  and  showed  that  the  mass  action 
laws  were  closely  related  to  a  least-action  principle  with 
respect  to  that  function.  That  is,  he  shov/ed  that  under 
certain  circumstances  minimizing  the  free  energy  function 
is  equivalent  to  satisfying  the  mass  action  laws. 

Until  recently,  this  relationship  between  the  mass 
action  laws  and  the  Gibbs  free  energy  function  was  ordi¬ 
narily  used  only  for  theoretical  purposes,  while  the 
actual  computation  of  equilibrium  compositions  normally 
employed  the  mass  action  laws.  It  has  recently  been 
shown  [6-8]  that  chemical  equilibrium  compositions  could 
be  determined  numerically  by  minimizing  the  Gibbs  free 

^See,  for  example.  Refs.  2-5. 
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energy  function.  Tb  .  technique  appears  to  offer  extremely 
powerful  and  far  reaching  tools  for  the  analysis  of  complex 
chemical  systems.  It  has  been  employed  in  the  analysis  of 
extremely  complex  models  of  physiological  systems  which 
were  formerly  far  beyond  the  reach  of  any  quantitative 
analysis  [9-12]. 

Due  to  the  fact  that  it  is  now  possible  to  analyze 
the  equilibrium  behavior  of  very  complex  chemical  systems 
by  means  of  automatic  digital  computation,  certain  ques¬ 
tions  which  appear  to  have  been  left  dormant  since  the 
time  of  Gibbs  become  important.  We  have  in  mind  such 
question.^  as:  Do  the  equations  given  by  the  mass  balance 
laws  and  mass  action  laws  possess  a  solution?  Is  this 
solution  unique?  If  it  is  not  unique,  what  is  the  nature 
of  the  set  of  all  solutions?  Does  the  Gibbs  free  energy 
function  possess  a  minimum?  Is  this  minimum  attainable 
by  a  unique  composition?  If  not,  what  is  the  nature  of 
the  set  of  all  minimizing  compositions?  And  finally, 
what  are  the  relations  between  the  solutions  to  the  mass 
action  laws  and  the  minima  of  the  Gibbs  free  energy  func¬ 
tion? 

These  questions  now  become  important  because,  first . 
the  extremely  large  chemical  systems  which  are  now 
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susceptible  to  analysis  are  more  apt  to  exhibit  pathologies, 
and  are  apt  to  exhibit  more  complex  pathologies,  than  the 
simpler  systems  with  which  we  have  been  able  to  deal  in  the 
past;  and  secondly,  the  use  of  automatic  digital  computa¬ 
tion  requires  a  much  more  precise  knowledge  of  the  nature 
of  the  pathologies  which  can  occur. 

In  Part  I  (Sections  1-5)  we  will  motivate  formulations 
of  these  problems.  In  Part  II  (Sections  6-12)  we  will 
attempt  to  answer  these  questions,  using  mathematical  tools 
which  were  unavailable  to  Gibbs. 

An  attempt  has  been  made  to  write  Part  I  so  that  it 
will  be  readable  by  the  mathematically  literate  chemist, 
physiologist,  etc.  Part  I  also  serves  as  a  chemical  in¬ 
troduction  to  the  mathematician.  Unfortunately,  the  chemist 
will  find  the  chemical  exposition  of  Part  I  quite  elementary 
and  belaboring  of  details.  The  mathematician  will  have  an 
analogous  complaint  as  to  the  mathematical  exposition  of 
Part  I. 

Part  II  will  require  some  familiarity  with  mathematical 
analysis,  linear  algebra,  and  the  theory  of  convex  functions 
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PART  I 


1.  THE  CHEMICAL.  SYSTEM 

We  will  consider  chemical  systems  composed  of  a  finite 
number  of  homogeneous  phases.  3y  a  homogeneous  phase  we 
mean  one  which  is  homogeneous  in  chemical  composition, 
pressure,  and  temperature.  The  question  of  under  what 
circumstances  a  phase  is  to  be  regarded  as  homogeneous  is 
one  which  is  not  considered  here.  It  can  only  be  answered 
as  part  of  the  process  by  which  we  formulate  a  model  of  an 
actual  chemical  system. 

For  example,  a  chemical  system  consisting  of  a  liquid 
solution  in  a  vessel  over  which  there  is  a  vapor,  might 
be  regarded  as  having  tv/o  phases--a  gas  phase  and  a  liquid 
phase.  A  chemical  system  consisting  of  two  solutions 
separated  by  a  semipermeable  membrane  might  be  regarded 
as  having  two  liquid  phases. 

It  should  be  emphasized  that  all  that  is  required  of 
a  phase  is  that  it  be  homogeneous.  It  need  not ,  for  example, 
occupy  contiguous  portions  of  space.  Thus,  in  Ref.  12,  all 
the  interiors  of  the  red  cells  of  the  body  are  usefully 
regarded  as  forming  a  single  phase.  (Note  that  this  phase 
is  not  really  homogeneous--but  no  phase  is  ever  really 
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honogeneous .  The  only  question  is:  Does  the  assumption 
of  homogeneity  hold  to  a  sufficient  degree  for  the  analyt¬ 
ical  purposes  which  the  designer  of  a  model  has  in  mind?) 

In  what  follows,  the  adjective  "homogeneous"  is  generally 
omitted  in  referring  to  homogeneous  phases.  Note  that  the 
possibility  of  a  chemical  system  consisting  of  exactly  one 
phase  has  not  been  precluded. 

It  shall  be  assumed  that  each  phase  is  at  some  specified 
pressure  and  temperature. 

Let  a  chemical  system  have  K  phases,  P,  , P2  , .  .  .  , ,  •  •  .  > 

Let  the  chemical  species  which  can  occur  in  the  phase, 

be  ) • • • ,  where  is  the  number  of  species 

k 


in  phase  k. 

This  last  statement  requires  some  elaboration.  First, 
by  a  species  is  meant  a  molecular  species.  Examples  of 
species  are  H2O,  NaCl,  Na  .  Note  that  when  a  species  is 
specified,  a  molecular  or  ionic  structure  is  implied  as 
well  as  its  composition  in  terms  of  the  number  of  atoms 
of  each  element  composing  a  molecule  of  the  species.  Thus, 
two  distinct  species  might  be  identical  when  looked  at  in 
terms  of  their  molecular  formulas.  Also,  for  reasons  cf 
technical  mathematical  convenience,  a  species  which  can 
occur  in,  say,  two  different  phases  will  be  regarded  as 
representing  two  different  species.  Thus,  if  H2O  could 
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occur  in  both  phase  and  phase  P2  ,  we  would  regard  H2O 
in  phase  P^,  and  H2O  in  phase  representing  two 

different  species.  Thus,  if  P^  is  a  liquid  phase  and  P2 
a  gaseous  phase,  we  might  refer  to  H20-liquid,  and  to  H2O- 
vapor.  When  a  molecule  moves  from  one  phase  to  another 
(for  example,  by  evaporation,  condensation,  or  migration 
across  a  semi- permeable  membrane)  we  will  regard  a  type  of 
chemical  reaction  as  having  occurred--a  reaction  in  which 
one  species  was  transformed  into  another  species. 

Secondly,  by  saying  that  a  species  can  occur  in  a 
phase  Pj^,  it  is  not  necessarily  meant  that  at  equilibrium 
any  positive  amount  of  will  be  found  in  Pj^.  It  is 
merely  that  provision  is  to  be  made  for  the  possibility 
that  Aj^^  might  occur  in  P^.  The  question  of  what  species 
will  be  provided  for  in  what  phases  is  a  question  which 
can  only  be  answered  by  judgment  and  experience,  and  by 
some  knowledge  of  the  characteristics  of  the  system  under 
study.  Theoretically,  one  might  provide  for  the  possible 
presence  of  each  of  the  billions  of  chemical  species  in 
each  phase.  If  this  were  done,  and  if  the  resulting 
numerical  problem  could  be  solved,  almost  all  species 


would  occur  in  negligibly  small--less  than  trace--amounts . 
Thus,  the  problem  of  prescribing  which  species  are  to  be 
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provided  for  in  which  compartments  is  a  problem  of  pre¬ 
dicting  which  species  will  occur  in  more  than  negligible 
amounts.  Of  course,  in  case  of  doubt,  one  can  provide  for 
a  few  species  for  which  one  is  not  certain. 

It  should  be  observed  that  if  a  phase  is  a  chemical 
solution,  then  the  solvent  itself  is  one  of  the  species 
occurring  in  P^. 

Let  be  the  number  of  moles  of  present  in  P^^. 
That  is,  the  number  of  molecules  of  present  is  (by 
definition  of  a  mole)  equal  to  multiplied  by  Avagadro's 
number  (approximately  10^^) . 

Let  N|^  be  a  vector  whose  components  are  >  •  •  •  > * 

k 

Let  N  be  the  vector  whose  components  are 

(Nil  ’^12  ’  •  *  *  ’^lli ’^^21  ’^22  ’  "  '  ’^212  ’  *  "  *’\l  ’^K2  ’  *  ‘ 


Thus,  N^  is  a  vector  with  components  which  specifies 
the  composition  of  the  phase  P^,  and  N  is  a  vector  with 


K 

^  “  I 

k=l 


components  which  specifies  the  composition  of  the  entire 


system. 


-9- 


Since  the  temperatures  and  pressures  of  the  phases 
are  assumed  specified,  the  vector  N  may  be  regarded  as 
a  state  vector.  That  is,  N  completely  describes  the  state 
of  the  system. 

Let  be  the  total  number  of  moles  in  when  the 

system  is  in  state  N;  that  is,  let 

(1.1)  a|^(N)  =  ^  .  for  1  i  k  s  K  . 

i=l 


Also ,  let 


for  l^^k^N,  l^i^I,  . 

k 


That  is,  is  the  concentration  of  when  the 

system  is  in  state  N.  Concentration  can  be  measured  in 
a  variety  of  scales.  Our  definition  of  n^^  amounts  to 
choosing  the  "mole  fraction"  scale. 


2.  THE  MASS  BALANCE  LAWS 

Following  Ref.  10,  we  will  express  the  mass  balance 
laws  as  linear  conditions  on  N.  Let  , . . .  ,B^  , .  . .  be 
a  set  of  fundamental  building  blocks .  such  that  each  species 
Ap.^  is  composed  of  these  building  blocks. 
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For  example,  , . .  .  ,B ,  , . . .  might  be  taken  to  be  the 

atomic  elements.  However,  it  is  often  better  to  choose  more 

complex  structures  as  building  blocks.^  Let  each  molecule 

of  A,  ,  contain  S,  .  ,  units  of  B,.  Thus,  the  vector 
ki  k,i,^  f  ’ 

(S,  .  S,  .  -)  is  essentially  the  molecular  formula 

for  Note  there  is  nothing  to  prevent  two  different 

species  from  having  identical  molecular  formulas.  Let 
be  the  total  number  of  units  of  B^,  in  the  system.  We 
than  have  the  conditions: 

(2.1)  ^  ^  ^k,i,^  \i  "  *  for  1  ^  -  L  . 

k=l  i=l 

The  description  of  those  aspects  of  the  system  con¬ 
cerned  with  such  questions  as  which  species  are  permeable 
to  interphase  boundaries ,  can  be  implemented  by  an  ap¬ 
propriate  extension  of  the  set  of  fundamental  building 
blocks.  This  process  has  been  more  completely  described 
elsewhere  llOj;  we  will  content  ourselves  here  with  a  simple, 
brief  example.  Consider  a  system  composed  off  two  phases, 
separated  by  a  semipemieable  membrane.  Suppose  each  phase 
contains  the  following  species:  H2O,  h"^,  OH  ,  a  protein, 

R,  and  an  ionized  form  of  that  protein,  RH^.  Suppose  the 

^See  Refs.  10  and  11. 
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membrane  is  permeable  to  all  species  except  R  and  RH  . 

This  situation  can  be  expressed  by  taking  as  building 

blocks ,  H',  OH  ,  and  two  forms  of  R:  R^,  representing 

R  in  phase  1;  and  R2,  representing  R  in  phase  2.  Thus, 

the  following  tableau  gives  the  values  of  the  S,  .  ^ . 

The  columns  represent  building  blocks;  the  rows  represent 

species;  the  body  of  the  table  gives  the  S,  .  . 

tc ,  i ,  I 


Table 

I 

h'*' 

oh" 

h 

«2 

"1" 

Ai  ^  =  H  in  phase  1 

1 

0 

0 

0 

Ai  2  =  OH  in  phase  1 

0 

1 

0 

0 

Ai  3  =  H2O  in  phase  1 

1 

1 

0 

0 

A^  ^  =  R  in  phase  1 

0 

0 

1 

0 

A^  3  =  RH*^  in  phase  1 

1 

0 

1 

0 

A2  ^  =  H'  in  phase  2 

1 

0 

0 

0 

A2  2  ~  OH  in  phase  2 

0 

1 

0 

0 

A2  3  =  H2O  in  phase  2 

1 

1 

0 

0 

A2  ^  =  R  in  phase  2 

0 

0 

0 

1 

A2  3  =  RH*^  in  phase  2 

1 

0 

0 

1 

We  now  leave  the  example  and  return  to  our  discussion 
of  general  chemical  systems. 
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A  chemical  system  may  have  one  or  more  phases  on 

which  it  is  desirable  to  impose  a  condition  of  charge 

neutrality.  To  impose  such  a  condition  on,  say,  phase 

P,  ,  let  S.  ,  T  be  the  charge  per  mole  of  A.,  .  (S.  ,  ^11 

k’  i,k,L+l  b  f  ik  i,k,L+l 

will  be  zero  unless  is  an  ion.)  Let  =  0 

I  /  0 

for  k  f  k.  Let  =  0.  Then,  the  charge  neutrality 

condition  for  has  the  same  form  as  the  equations  of 
(2.1)  and  can  be  included  among  them. 

Finally,  from  the  definition  of  it  is  clear  that 

we  must  have 


(2.2)  \i  ^  0  1  ^  k  ^  K  ,  1  ^  i  ^  . 

We  are  hence  led  to  view  the  constraints  on  the  system 
as  a  set  of  simultaneous  linear  equations  of  the  form 

(2.1)  and  a  set  of  simultaneous  inequalities  of  the  form 

(2.2) . 

3.  REACTION  VECTORS  AND  THE  STOICHIOMETRIC  CONDITIONS 
We  consider  chemical  reactions  involving  the  species 
of  a  system,  in  equilibrium  chemistry,  a  reaction  is 
described  by  writing  two  formal  sums  separated  by  a  double 
arrow.  Each  of  these  formal  sums  consists  of  one  or  more 
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terms.  Each  of  these  terms  consists  formally  of  a  positive 
real  number,  called  a  stoichiometric  coefficient,  multi¬ 
plied  by  the  S3mibol  for  some  chemical  species  of  the  system. 
The  species  occurring  in  the  left-hand  formal  sum  are  called 
reactants;  the  species  occurring  in  the  right-hand  formal 
sum  are  called  products.  The  two  formal  sums  separated  by 
a  double  arrow  form  the  stoichiometric  equation  of  the 


reaction.  Note  that  mathematically  speaking,  the  stoichi¬ 
ometric  equation  is  not  an  equation  at  all. 

Given  any  such  stoichiometric  equation,  we  shall  define 
a  vector  9,  with  components  0^^^  (l^k^K,  l^is^  I^) , 
as  follows.  To  define  9,  consider  the  reaction  as  written 
out.  If  does  not  occur  in  the  reaction,  either  as  a 
product  or  a.s  a  reactant,  let  0^^^  "0.  If  occurs  as 
a  reactant,  but  not  as  a  product,  let  9^^^^^  be  its  stoichi¬ 
ometric  coefficient.  If  occurs  as  a  product,  but  not 
as  a  reactant,  let  0^^^^^  be  the  negative  of  its  stoichiometric 
coefficient.  Finally,  if  occurs  both  as  a  product  and 
as  a  reactant,  let  9^^^  be  its  coefficient  as  a  reactant, 
minus  its  coefficient  as  a  product. 

To  understand  the  definition  of  8,  treat  the.  stoichi¬ 
ometric  equation  formitlly  as  a  mathematical  equation, 
bringing  all  terms  to  the  left-hand  side  (changing  the 
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signs  in  the  process)  and  collecting  terms.  The  vector 
r  has  been  defined  so  that  this  manipulation  yields  the 
’'equation" 

K 

I  I  — 

k=l  i=l 

We  shall  call  any  vector  9  so  obtainable  from  a  re¬ 
action  a  reaction  vector.  Observe  that  different  re¬ 
actions  may  have  the  same  reaction  vector. 

We  have,  as  a  condition  for  a  vector  6  to  be  a  re¬ 
action  vector,  that 

K 

I  I  °  •  for  1  -  .  *  L  , 

k»*l  i”l 

where  the  ^  those  of  equations  (2.1).  The  con¬ 

dition  is  sufficient  as  well  as  necessary;  we  shall  re¬ 
gard  any  vector  6  that  satisfies  (3.2)  as  &  reaction  vector. 

Note  that  (2.1)  end  (3.2)  are  closely  related^  They 
differ,  however,  in  two  respects:  1)  in  that  (2.1)  is  a 
condition  on  a  vector  N  in  order  that  it  be  a  composition 
vector,  while  (3.2)  is  a  condition  on  a  vector  0  in  order 
that  it  be  a  reaction  vector;  and  2)  in  that  the  right-hand 


-15- 


sides  of  (3.2)  are  all  zero,  while  those  of  (2.1)  are  not. 

Unfortunately,  the  chemical  literature  refers  to  both 
sets  of  conditions  as  mass  balance  laws.  We  shall  reserve 
the.  term  ''mass  balance  laws"  for  (2.1)  and  call  (3.2)  the 
stoichiometric  conditions. 


.  THE  MASS  ACTION  LAWS 


We  shall  confine  our  attention  to  systems  which  ex¬ 
hibit  the  simplest  form  of  the  mass  action  laws.  Such 
systems  are  called  "ideal"  in  Ref.  13. 

In  our  terminology,  the  mass  action  laws  can  be  stated 
as  follows:  For  any  reaction  having  a  reaction  vector  9, 
there  is  an  equilibrium  constant  k(e),  such  that,  for  any 
composition  vector  N  representing  an  equilibrium  state  of 
the  system,  we  have 


K  \  ‘ 

(4.1)  n  n  (tv.(N))  =  k(5)  . 

k=l  i=l 


We  shall  find  it  useful  to  perform  certain  manipula- 

I 

tions  on  (4.1).  Let  us  first  take  the  logarithm  of  each 
side  to  get 


The  results  of  this  paper  also  apply  to  certain  types 
of  non- ideal  systems.  However,  we  prefer  to  simplify  the 
exposition  by  expressing  our  results  in  terms  of  ideal  systems 
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I  I  ®ki  ' 

k=l  i=l 

Equation  4.2  is  thus  satisfied  for  all  vectors  8  which 
are  reaction  vectors;  that  is,  for  all  vectors  0  v^ich 
satisfy  the  stoichiometric  condition  (3.2).  It  should  be 
noted  that  the  form  of  (3.2)  implies  that  any  scalar 
multiple  of  a  reaction  vector  is  a  reaction  vector,  and 
that  the  sum  of  any  two  reaction  vectors  is  a  reaction 
vector.  Furthermore,  the  form  of  equation  (4.2)  shows 
that  if  0  and  cp  are  reaction  vectors,  then 


log  k(8  +  cp)  =  log  k(8)  +  log  k(cp)  , 


and  that  if  r  is  a  scalar,  then 

log  k(r0)  =  r  log  k(e)  . 

These  facts  imply^  the  existence  of  constants 
that 


such 


(4.3)  log  k(e) 


I 


k=l 


k 

I  Si  'ki  • 

i=l 


Be  <,ause  log  k(0)  is  then  a  linear  function  of  0. 
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The  substitution  of  the  value  of  log  k(8)  from  (4.3) 
into  (4.2)  allows  us  to  restate  the  mass  action  laws  as 
follows:  There  are  constants  (l^k^K,  1^:  i^  I^) 
such  that  if  6  is  any  reaction  vector,  and  if  N  is  any 
equilibrium  composition,  then 

K  ^k 

I  I  \i  ("ki  +  "ki(«»  =  0  • 

k=l  i=l 

This  form  of  the  mass  action  law,  which  is  quite  equiv¬ 
alent  to  the  standard  forms,  is  more  convenient  for  our 
purpose. 

Another  form  is  often  quite  useful  for  computational 
purposes.  Although  there  is  no  need  for  it  in  the  present 
paper,  it  is  developed  here  for  the  sake  of  completeness; 
the  reader  may  omit  the  rest  of  this  section  without 
destroying  any  continuity  with  what  follows  it. 

Let  us  define 

(*•5)  yki  =  '=ki  +  iVi . 


Let  us  think  of  y  as  an  I-dimensional  vector  where 
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K 


k=l 


If  we  do  this,  then  the  condition  (4.4)  states  that  y  is 
orthogonal  to  all  vectors  9  which  satisfy  the  stoichio¬ 
metric  equations  (3.2).  If  we  think  of  ths  quantities 
S,  .  ,  as  forming  L  vectors  S, , . . .  ,S  , . . .  ,5,.  ,  each  having 
dimension  I,  then  (3.2)  becomes:  9  is  orthogonal  to  each 
S.  Thus,  the  mass  action  law  may  be  restated  as  follows: 
y  is  orthogonal  to  every  vector  which  is  orthogonal  to  each 
S.  But,  it  follows  from  elementary  linear  algebra  that 
the  condition  is  satisfied  if  and  only  if  y  is  representable 
in  the  form 

(4.6)  y  =  ^  TT^ 
or,  equivalently, 

L 

(4.7)  +  log  ^  . 

t-1 

Thus,  the  satisfaction  of  the  mass  action  laws  is 
equivalent  to  the  existence  of  quantities 
which  satisfy  (4,7)* 


-19- 


5.  THE  GIBBS  FREE  ENERGY  FUNCTION 

We  will  have  frequent  occasion  to  refer  to  the  free 
energy  of  a  system.  Since  we  are  dealing  with  ideal 
systems  and  since  pressure  and  temperature  are  constant 
for  each  phase,  the  free  energy  of  a  system  is  a  function 
F(N)  of  its  composition  N.  It  may  be  written  as 
follows  [13,14] 

K 

(5.1)  F(N)  =  ^  ^  N^.  (C^.  +  log  \i(N)) 

k=l  i=l 

where  (a  constant)  is  called  the  free  energy  parameter 
of 

The  least-action  principle  for  the  free  energy  states 
that  the  system  is  in  equilibrium  if  and  only  if  its  free 
energy  is  a  minimum,  subject  to  the  mass  balance  con¬ 
straints  (2.1)  and  to  the  constraint  (2.2);  that  is  to 
say,  a  composition  vector  represents  an  equilibrium 
state  if  and  only  if  F(N^)  ^  F(N)  holds  for  all  states  N 
which  satisfy  the  conditions  (2.1)  and  (2.2)  and  which 
are  sufficiently  close  to  N^. 

Note  tha^  the  least-action  principle,  as  stated  here, 
requires  only  that  F  have  a  local  minimum  at  N^.  It  will. 
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however,  turn  out  that  every  local  minimum  of  F  is  really 
a  global  minimum,  and  hence  that  the  underlined  phrase 
above  can  be  eliminated. 
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PART  II 


6 .  NOTATION 

Let  E  be  the  real  line.  For  any  integer  J  let  E’^  be 

the  real  Euclidean  space  of  dimension  J.  For  a  point 

N  e  E*^  we  will  say  that  N  -  0  when  all  the  components  of 

N  are  non-negative.  We  will  say  that  N  >  0  when  all  the 

components  of  N  are  positive. 

Let  K  be  a  positive  integer  representing  the  number 

of  phases.  Let  , .  .  . , . . . , be  positive  integers, 

th 

being  the  number  of  species  in  the  k  phase.  Let 


I  = 


We  will  find  it  convenient  to  designate  the  components  of 
a  point  N  €  E^  by  means  of  a  doubly- indexed  system  of  sub 
scripts.  Thus: 


N  = 


) 


KI 


K 


For  each  k  (1  ^  k  ^  K),  let  denote  the  I^-dimensional 
vector  whose  components  are  (1  ^  i  ^  ^k^  ’  thus, 
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)  . 


th 

then  represents  the  composition  of  the  k  phase.  We 
can  then  symbolically  write 


N  )  •  •  •  >  •  •  •  * 

7 .  QUASI-DEPENDENCE 

DEFINITION  7.1.  Two  points  M,  N  ?  are  quasi¬ 
dependent  if  for  each  phase,  k,  the  vectors  and  N,^ 
are  linearly  dependent;  that  is,  if  for  each  k  there 

exist  scalars  a,  and  b,  ,  not  both  zero,  such  that 

k  k 

+  bj^N^  =0.  If  M  and  N  are  quasi-dependent,  we 
will  write  M  ~  N. 

The  following  properties  of  quasi-dependence  are 
readily  apparent  from  the  definition. 

THEOREM  7.2  Let  M,  N,  N'  e  E^. 

A.  M  ~  M. 

B.  If  M  ~  N,  then  N  M. 

C.  If  M  and  N  are  quasi-dependent,  then  any 

two  linear  combinations  of  M  and  N  are  quasi¬ 
dependent  . 

D.  If  M  ~  N  and  M  ~  N'  ,  then  M  N  N'  . 
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We  note  that  is  not  a  transitive  relation: 

N  ~  M  and  M  N'  do  not  necessarily  entail  N  ~  N'-- 
for  example,  if  M  =  0. 

THEOREM  7.3.  Let  M,  N  e  E^.  A  necessary  and 
sufficient  condition  for  M  ~  N  is  that  for  every  phase 
k  and  all  i  and  j  with  1  ^  i,  j  ^  we  have 


(7.1)  M^j  . 

Proof.  First,  to  prove  sufficiency,  assume  that 
(7.1)  holds  and  let  k  be  any  phase.  If  =  0,  then 
and  are  clearly  linearly  dependent.  Otherwise,  there 
is  a  j  with  0.  With  j  regarded  as  fixed  and  i  as 

variable,  (7.1)  then  says 


N,  =  0  . 


But  since  4  0,  this  shows  that  and  are  linearly 
dependent . 

Second,  to  prove  necessity,  let  M^  and  be  linearly 
dependent.  There  are  constants  aj^  and  b^,  not  both  zero, 
with  =0.  We  may  assume  without  loss  of 

generality  that  a^^  ^  0.  We  then  have  for  all  i 
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M,  .  =  -  —  N,  .  . 
Ki  a,  ki 
k 


Hence , 

«ki 


—  N,  .  N,  . 
a^  ki  kj 


N,  .  N,  .  •■=  M,  .  N,  . 
a,  kj  ki  kj  ki 


Q.E.D. 


DEFINITION  7.4.  Let  M,  N  E  .  We  will  say 

that  M  is  quasi-dependent  on  N  it,  for  each  phase  k, 

the  vector  is  a  scalar  multiple  of  N^.  For  any 

I 

N  e  I,  let  /'(N)  denote  the  set  of  all  M  E  which 
are  quasi-dependent  on  N. 

We  thus  see  that: 

THEOREM  7.5.  A(N)  is  the  set  of  all  points  in 
E^  which  can  be  written  in  the  form 


(7.2)  (X^N3^,...,X^N^,...,X^N^)  , 

where  the  X^^  are  scalars.  A(N)  is  a  linear  subspace 
of  E".  The  dimension  of  A(N)  is  the  number  of  phases 
k  for  which  N^  0.  If  M,  M*  c  A(N),  then  M  ~  M'. 
Finally,  M  e  A(N)  if  and  only  if  M  ~  N,  and  M^^  =  0 
whenever  N^  =  0. 

Proof.  Clear. 

The  following  lemma  will  be  used  in  a  later  proof. 
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LEMMA  7.6.  Let  N  s  0,  M  e  A(N) ,  M  ^  0,  and 

M  7^  N.  Then  there  exists  M*  e  A(N)  such  that  M*  ^  0 

and  M'  ^  N,  and  such  that  N  lies  on  the  interior  of 

the  line  segment  joining  M  to  H' . 

Proof.  Let  X,  ,...,Xt,  be  the  scalars  in  the 

-  1  k  K 

representation  (7.2)  of  M  e  A(N).  Since  M  ^  N,  not  all 
the  X^  are  1.  Hence,  the  number  c  =  maXj^|l-X^|  is  posi¬ 
tive.  Define  X^  =  l-(X^-l)/c,  and  define  M^^^  =  ^k^k’ 
each  k.  Then  M'  e  A(N) .  Now,  (Xj^-l)/c  never  exceeds  1, 
and  is  nonzero  for  at  least  one  k.  It  follows  that  M'  ^  0 
and  M'  N.  Moreover,  we  have 

N  =  M  +  ^  M'  .  Q.E.D. 

1+c  1+c  ^ 

Note.  To  understand  the  relevance  of  quasi-dependence, 
let  us  consider  its  meaning  as  applied  to  two  points,  M, 

N  e  E^  with  M  >  0  and  N  >  0.  Then,  by  Theorem  7.5,  M  is 
quasi-dependent  on  N  if  and  only  if  M  and  N  are  quasi- 
dependent.  Furthermore,  in  the  notation  of  Sec.  1,  we 
have  (by  Theorem  7.3)  that  M  ~  N  if  and  only  if  n(M)  =  n(N) 
that  is,  M  ~  N  if  and  only  if  M  and  N  are  composition 
vectors  yielding  identical  concentrations. 
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8.  PROPERTIES  OF  THE  FREE  ENERGY  FUNCTION 

The  free  energy  function  as  given  by  (5.1)  is  de¬ 
fined  only  for  those  N  e  for  which  N  >  0.  This  is  so 
because  log  t  is  undefined  when  t  ^  0.  However,  the 
definition  can  be  extended  to  all  N  e  for  which  N  ^  0. 
To  do  this,  first  let  us  repeat  the  definition  (1.1): 

DEFINITION  8.1.  Let  1  ^  k  ^  K  and  let  N  e  E^. 

Then 

=  I  \i  • 

i=l 

If  we  expand  (5.1),  taking  note  of  (1.2),  and 
simplify  the  result,  we  arrive  at  the  extended  defini¬ 
tion  of  F(N) ; 

DEFINITION  8.2. 

K  '•k  K 

F(N)  =  C-N  +  y  ^  log  ^  o^(N)  log  o^(N)  , 

k=l  i=l  k=l 

where  0  log  0  is  taken  to  be  zero,  and  where 

Z^ik«ik  • 
i,k 
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This  then  defines  a  function  F(N)  for  all  N  0, 

because  if  N  ^  0,  then  0,  for  all  k.  Also, 

since  lim,  t  log  t  =  0,  we  have: 
t-0 

THEOREM  8.3.  F(N)  is  defined  and  continuous  for 


all  N  ^  0. 

We  will  wish  to  consider  directional  derivatives  of  F. 
Clearly,  it  does  not  make  sense  to  speak  of  the  derivative 
of  F  at  N  in  the  direction  unless  N  +  t  is  in  the 
domain  of  F  for  all  small  positive  t.  This  motivates  the 
following  definition: 

DEFINITION  8.4.  Given  -  ,  N  €  E^,  we  will  say 
that  ®  is  admissible  atNifN+t^'^O  for  all  suf¬ 
ficiently  small,  positive,  real  t. 

The  following  facts  about  admissibility  will  be 
useful : 


LEMMA  8.5.  Let  9,  N  e  E  .  A  necessary  and  suf¬ 
ficient  condition  for  0  to  be  admissible  at  N  is  that 


N  ^  0  and  that  .  -  0  whenever  N,  .  =  0. 

kj  kj 

Proof.  Clear. 


LEMMA  8.6.  Let  8  be  admissible  at  N.  For  each 

k  (1  ^  k  ^  K)  let  be  the  set  of  all  j  for  which 

either  N,  .  0  or  9  ^  ^  0.  Let  J  be  the  set  of  all 

kj  ^  kj 

k,  for  which  either  ^  0  or  9^  0;  then: 
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8 .6A:  Jj^  is  empty  for  all  k.  for  which  k  ^  J. 

8 . 6B:  N,  ,  +  t  8,  ,  >  0  for  all  k,  j,  with  j  e  J  , 
kj  kj  ic 

and  for  all  sufficiently  small,  positive  t. 

Proof  of  8.6A.  Let  k  J.  We  have,  by  the  definition 
of  J, 

1)  a^(N)  =  0 

2)  =  0. 

By  Lemma  8.5,  N  ^  0.  Hence,  using  (1): 

3)  =  0. 

Hence,  by  Lemma  8.5,  B.  a  Q;  but  then,  using  (2): 

4)  9^=0. 

But  (3)  and  (4)  imply  that  is  empty.  Q.E.D. 

Proof  of  8.6B.  If  >  0,  the  result  clearly  holds. 

Otherwise  (since  N  ^  0  by  Lemma  8.5),  N,  .  =  0,  and  by 

kJ 

Lemma  8.5,  9^^^  ^  0;  but  by  the  definition  of  Jj^,  v/e  see 

that  N,  .  and  9,  .  cannot  both  vanish.  Hence,  9,  .  >  0,  and 
k J  kj  KJ 

again  the  result  holds. 

LEMMA  8.7.  Let  9  be  admissible  at  N.  Let  the 


sets  J,  ,  and  the  set  J  be  defined  as  in  Lemma  8.6. 
k 

Let  2  =  N  +  t9.  Then,  for  all  sufficiently  small 


t  >  0,  we  have: 
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8.7A: 


F(Z)  -C.Z+  I  ^  Z^.  log  Z^.  -  I  a^(Z)  log  ,^(Z)  ; 

kcj 

and , 

8 .7B;  F(N  +  t9)  is  differentiable  with  respect  to 
t ;  in  fact 


±_ 

dt 


F(Z)  = 


E 

k-rJ  jeJ. 


kj 


log  Z 


kj 


keJ 


> 


and, 

8.7C; 


d_ 

dt 


F(Z)  = 


and , 

8 .7D;  F(N  +  t&)  is  twice-differentiable  with  respect 
to  t;  in  fact,^ 


^This  formula,  for  the  case  K  «=  1  and  N  >  0,  appeared 
in  Ref.  6. 
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Aj-so,  we  have, 


LIE;  If  for  some  k  and  j  „a  have  =  o,  o^(N)  >  0, 

®kj  iE  E(Z)  =  -  ».  Finally. 

If  for  all  k  and  j  for  which  N  =  0  and 

\W  >  0,  we  have  9^.  =  0  (i.e..  if  the  hypothesis 
of  8.7E  is  not  satisfied),  then 


^  dl  “ 


'A 


+  log 


keJ  jej 

N  >0 
kj 


^^(N) 


S£2o£.o£_8JA.  This  formula  follows  from  Definitio 
8.2  and  the  convention  concerning  0  log  0. 

If  we  differentiate  both  sides  of 


8.7A,  we  obtain 
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dF(Z) 

dt 


C-  +  ^  I  (1  +  los  \j)  -  I  (1  +  1°8 

keJ  kej 


This  differs  from  the  right  side  of  8.7B  by 


1  Z  ®kj  ■  E  • 

keJ  icJ,  keJ 

k 


But  this  is  seen  to  be  0,  using  Definition  8.1  and  the 
fact  that  9,  .  =  0  when  j  /  J,  . 

rCJ  R 

Proof  of  8.7C.  Follows  from  8.7B.  The  right-hand 
side  is  well  defined,  because  of  8.6B. 

Proof  of  8.7D.  If  we  differentiate  both  sides  of 
8 .7B  we  get: 


(1) 

'  dt 


z 

ke  J 


<■ 

V 


z 

k€j 


a^(Z) 


But  if  we  expand  the  inner  term  of  the  right-hand 


side  of  8.7D,  we  get 
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( right -hand  side  of  8.7D)  = 


=  (right-hand  side  of  (1)).  Q.E.D. 

Proof  of  8»7E  and  8.7F.  Both  of  these  propositions 
follow  from  8,7C,  as  can  be  seen  by  separating  the  sum¬ 
mation  of  8.7C  into  three  portions,  depending  on  the 
following  conditions  on  k  e  J  and  j  e 
Portion  1:  N,  .  >  0. 

Portion  2:  =  0  and  =  0. 

Portion  3:  =  0  and  >  0. 


Portion  3  will  be  vacuous  in  the  case  hypothesized 
by  8.7F  and  Portions  1  and  2  converge  to  the  two  sums  of 
8.7F.  If  Portion  3  is  not  vacuous,  it  converges  to  -■», 


Q.E.D. 
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DEFINITION  8.8.  Let  -  be  admissible  at  N.  De¬ 
fine  the  directional  derivative  F^(N)  to  be 


def  ii_  F(N  9-  t  6)  -  F(N) 
t-0+ 


in  the  event  that  this  limit  exists,  accepting  +“ 
and  -®  as  possible  limits. 

In  order  to  evaluate  Fg(N),  we  will  need  the  follow¬ 
ing  consequence  of  the  mean  value  theorem  of  elementary 
calculus . 

LEMMA  8.9.  Let  t^  >  0  and  let  g(t)  be  a  con¬ 
tinuous  real-valued  function  in  the  interval  0  ^  t  ^  t  . 

o 

Suppose  that  exists  in  0  <  t  <  t^  and  that 


a  =  lim, 


dR(t:) 

dt 


exists  (accepting  +“  and  as  v.alues  for  a);  then, 


exists  and  is  equal  to  a. 
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THEOREM  8.10.  Let  5  be  admissible  at  N.  Sup¬ 
pose  that  for  some  k  and  j  we  have  N,  .  =  0,  c,  (N)  >  0, 

kj  k 

and  >  0.  Then,  F^CN)  =  -“. 

Proof.  Lemmas  8.7E  and  8.9. 

THEOREM  8.11.  Let  9  be  admissible  at  N.  Sup¬ 
pose  that  for  all  k  and  j  for  which  =  0  but 
'k(N)  >  0,  we  have  9^,  =  0.  Then 


Fl(N)  =  ^ 
keJ 


keJ 


+  log 


Ai 


+  log 


Si,* 
—Hi 

^k(s)  J ' 


Prop f .  Lemmas  8.7F  and  8.9. 

The  reader  might  ask  why  we  have  not  dealt  directly 
with  the  partial  derivatives  of  F  in  order  to  calculate 
such  quantities  as  Fg(N),  or  indeed  in  order  to  examine 
directly  the  behavior  of  F.  The  answer  lies  in  an  in¬ 
spection  of  the  formula  of  (8.11),  which  is  nonlinear  in 
0,  as  long  as  any  of  the  components  of  N  vanish;  that  is, 
as  long  as  the  second  summation  of  Theorem  8.11  is  not 


vacuous . 
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The  techniques  for  studying  the  extreme  points  of  a 
function  of  several  variables  through  their  first-order 
partial  derivatives  depends  on  the  fact  that,  for  suf¬ 
ficiently  well  behaved  functions,  the  directional  deriva¬ 
tive  of  a  function  G(N)  in  the  direction  -  is  a  linear 
function  of  6,  and  hence, 


Gi(N)  = 

i 


aG(N) 

^N. 

1 


where  the  summation  index  ranges  over  the  components  of 
Many  of  the  complications  of  this  paper  arise  because 
this  relation  fails  for  F(N) . 

LEMMA  8.12,  Let  ?  and  --  be  admissible  at  N. 
Then  Fl (N)  =  -F^^(N) . 

Proof.  We  must  have  (because  of  Lemma  8.5)  =  0 

- — 

for  all  k  and  j  for  which  =  0.  Hence,  the  result 
follows  from  Theorem  8.11. 

THEOREM  8.13.  F  is  convex  on  its  domain. 

Proof.  This  follows  immediately  from  the  fact  that 
F  is  continuous  (Theorem  8,3),  and  from  Lemma  8.7D,  which 
guarantees  that  F  is  convex  in  the  interior  of  all  line 
segments  contained  in  its  domain  [15,. 


-36- 


LEMMA  8.14.  Let  0  ~  N,  let  N  e  E^,  and  let 
0  be  admissible  at  N.  Then,  F(N  +  t9)  is  a  linear 
function  of  t,  for  all  t  (positive,  zero,  or  negative) 
for  which  N  +  t9  -  0. 

Proof.  Let  g(t)  =  F(N  +  t0).  We  will  first  consider 
only  the  t  >  0  for  which  N  +  tB  ^0.  By  Lemma  8.7D, 
Theorem  7.2C,  and  Theorem  7.3,  we  have 


for  t  >  0.  Hence,  g(t)  is  linear  for  t  >  0.  Similarly, 
g(t)  is  linear  for  t  <  0,  in  the  event  that  -0  is  also 
admissible.  That  the  two  linear  functions  are  the  same 
is  due  to  the  fact  that  they  coincide  at  t  =  0  (Theorem 
8.3),  and  to  Lemma  8.12,  which  says  that 


V 


Q.E.D. 


Conversely; 

LEMMA  8.15.  Let  t,  t  be  real  numbers  different 
from  zero.  Let  M,  N  e  E^.  Suppose  that  M  ^  0,  N  ^  0, 
tM  +  "tN  ^  0  and  that  F(tM  +  tN)  =  tF(M)  +  tF(N) .  Then, 


M  ~  N. 
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Proof .  By  Theorem  7.2,  and  the  fact  that  F(tN)  = 
tF(N)  for  all  t  ^  0,  N  ^  0,  we  may  assume  without  loss 
of  generality--by  a  change  of  notation--that  t  >  0,  t  >  0, 
and  t  +  t  =  1.  Then  we  know,  by  the  convexity  of  F 
(Theorem  8.13),  that  F  is  linear  on  the  line  segment 
joining  M  to  N  [15j.  The  result  then  follows  from  Lemma 
8.7D  and  Theorem  7.3. 

9.  THE  MINIMUM  SET  OF  THE  FREE  ENERGY  FUNCTION 

The  results  of  Sec.  8  will  now  enable  us  to  specify 
the  nature  of  the  set  on  which  F  takes  its  minimum. 

The  mass  balance  laws  of  Sec.  2  have  the  effect  of 
confining  N  to  a  certain  linear  manifold  in  E^.  According¬ 
ly,  we  shall  be  interested  in  the  minimum  of  F  on  a  fixed 
linear  manifold,  H  of  E^,  subject  to  the  additional  con¬ 
straint  that  N  ^  0. 

DEFINITION  9.1.  Let  Q  be  the  set  of  all  N  e  E^ 
for  which  N  ^  0.  Let  H  be  a  fixed  linear  manifold 
(the  "mass  balance"  manifold)  of  E^.  Let  D  be  the 
minimum  set  of  F  on  H  H  Q.  That  is,  let  D  be  the 
set  of  all  N  <=:  H  n  Q  such  that  F(N)  ^  F(M)  for  all 
M  e  H  n  Q. 

Note.  We  have  now  introduced  all  the  concepts  neces¬ 
sary  to  define  a  chemical  equilibrium  problem.  We  might 
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thus  define  an  abstract  chemical  system  to  be: 

1)  An  integer,  K  >  0. 

2)  Integers  0. 

3)  A  linear  manifold  H  of  ,  such  that  H  "  Q  is  not 
empty,  where 


k 

4)  A  point  C  ?  E^. 

K  specifies  the  number  of  phases.  The  specify  the 
number  of  species  in  the  various  phases.  K  determines 
the  mass  balance  laws.  C  determines  the  parameters  of 
the  free  energy  function  (or  of  the  mass  action  laws) . 
THEOREM  9.2.  If  D  is  nonempty,  there  is  a 
point  N  <?:  D  such  that^  D  =  '(N)  H  "  Q. 

We  shall  obtain  Theorem  9.2  as  a  consequence  of 
several  lemmas. 

LEMMA  9.3.  D  is  convex. 

Proof .  This  follows  from  the  fact  that  F  is  convex 
(Theorem  8.13)  and  the  fact  f  at  the  minimum  set  of  a 
convex  function  is  convex.*  Q.E  D. 

^Recall  the  definition  of  ‘ (N)  in  Definition  7.4. 
*See  Ref.  15. 
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LEMMA  9.4.  F  is  constant  on  D. 


Proof.  This  follows  from  the  definition  of  D.  Q.E.D. 

LEMMA  9.5.  If  M  M'  '  D,  then  M  M' . 

Proof .  If  M  =  M'  ,  then  the  result  holds  because  M  M 
(Theorem  7.2A).  It  W  ^  M' ,  then  by  the  convexity  of  D 
(Lemma  9.3),  D  contains  the  interval  joining  M  to  M' . 

Hence,  since  F  is  constant  on  D  (Lemma  9.4),  F  is  constant 
on  the  interval,  and  hence  linear  on  the  interval.  Under 
these  circumstances,  we  must  have  (Lemma  8.15)  MM'.  Q.E.D. 

LEMMA  9.6.  Let  N  D;  then,  *(N)  ^  H  "  Q  -  D. 

Proof.  Choose  any  point  M  €  ‘ (N)  '  H  "  Q  and  choose 
M'  €  ‘ (N)  "  Q  according  to  Lemma  7.6,  so  that  N  is  interior 
to  the  line  segment  L  joining  M  and  M'.  We  have  M  ~  M' 
(Theorem  7.5);  hence,  M  ~  M'  -  M  (Theorem  7.2C),  hence, 

F  is  linear  on  L,  by  Theorem  8.14.  But,  L  c  H  r  Q,  be¬ 
cause  M,  N  €  H  and  M,  M'  €  Q.  If  F  were  not  constant  on 
L,  then  either  F(M)  <■  F(N)  or  F(M')  <  F(N)  ,  contradicting 
N  D.  Hence,  F(M)  =  F(N)  ,  and  M  e  D.  Q.E.D. 

LEMMA  9.7.  If  D  is  nonempty,  there  is  a  point, 


N  e  D,  such  that  if  for  any  k,  «  0,  then  for  all 
M  D  we  must  have  M,  =0. 


Prop f .  For  each  k  (1  ^  k  ^  K),  define  a  point  m''  ^  as 
follows:  If  M  D  can  be  found  with  M^  0 ,  take  M^^^  =  M; 
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if  not,  let  be  an  arbifary  element  of  D.  Then  the 

point 


N 


M 


(k) 


has  the  required  properties,  by  convexity  (Lemma  9.3)  and 
the  fact  that  D  c  Q  by  definition  of  D. 

Proof  of  Theorem  9.2.  Let  N  be  the  point  given  by 
Lemma  9.7.  By  Lemma  9.6,  .\(N)  ^  H  n  Q  c  D.  Let  M  e  D. 

By  Lemma  9.5,  M  ~  N.  By  the  way  in  which  N  was  defined, 
for  no  k  can  Nj^  be  zero  unless  is  also  zero.  Hence, 
by  Theorem  7.5,  M  e  A(N) .  We  have  proved  that  D  c  A(N) . 
But,  by  the  definition  of  D,  D  c  H  n  Q.  Hence, 

D  c  A(N)  H  n  q.  Thus,  D  =  A(N)  n  H  n  Q,  Q.E.D. 

THEOREM  9.8.  Let  N  e  D.  Suppose  that  for  some 
k  and  j,  we  have  =  0,  but  ^  0.  Then,  for  all 
M  e  H  "  Q,  we  must  have  **  0. 

Prop f .  Suppose  that  for  some  M  e  H  H  Q,  we  have 

>  0.  Since  N  ^  0  and  ^  0,  we  have  >  0.  Let 

S  =  M  -  N.  Then,  N  +  t6  e  H  H  Q  for  all  real  t  for  which 

0  ^  t  ^  1.  Therefore,  F(N  +  tS)  ^  F(N)  for  0  ^  t  ^  1. 

Furthermore,  >  0.  But  then  by 

Theorem  8,10,  F(N  +  t9)  is  a  strictly  decreasing  function 
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of  t  for  all  sufficiently  small  positive  t — this  con¬ 
tradicting  the  assumption  that  N  D.  Q.E.D. 

Note.  The  import  of  Theorem  9.8  is  roughly  that  every 
species  (j)  whose  presence  has  been  allowed  for,  must  be 
present  in  some  positive  amount  in  any  equilibrium  com¬ 
position,  unless  the  entire  phase  (k)  containing  that 
species  is  not  present.  More  precisely,  if  no  amount  of 
j  is  present  at  equilibrium,  then  either  N,  ,  must  be  zero 
by  virtue  of  the  mass  balance  laws  alone,  or  no  amount  of 
any  j '  in  the  same  phase  as  j  can  be  present  for  that 
equilibrium  composition. 

Of  course,  nothing  in  this  theorem  forbids  N^^  being 
0  for  some  equilibrium  compositions  and  non-zero  for 
other  equilibrium  compositions. 

Note .  Theorems  9.2  and  9.8  provide  a  complete  de¬ 
scription  of  the  possible  equilibrium  sets,  D.  As  a 
matter  of  fact,  we  will  prove  below  (Theorem  11.1)  that 
any  nonempty  set  D  satisfying  the  conditions  of  these 
theorems  is  capable  of  being  the  equilibrium  set  for  some 
chemical  system. 

It  follows  from  the  fact  that  F  is  a  convex  function, 
that  D  is  a  convex  set.  But,  Theorem  9.2  asserts  much 
more  than  the  convexity  of  D.  In  fact.  Theorem  9.2 
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represents  D  as  the  intersection  of  a  linear  manifold, 
(•(N)  ''  H),  with  Q. 

THEOREM  9.9.  If  H  n  Q  is  bounded,  then  D  is 
nonempty. 

Proof.  Since  H  Q  is  closed  and  F  is  continuous 
(Theorem  8.3),  F  must  attain  a  minimum  value  on  H  H  Q. 

Q.E.D. 

Note.  There  are  no  results  in  this  paper  which  deal 
with  the  question  of  under  what  conditions  D  is  nonempty 
when  H  0  Q  is  unbounded.  Although  the  question  appears 
to  be  interesting  mathematically,  it  appears  to  be  of 
little  chemical  interest,  because  in  actual  applications 
H  n  Q  is  always  bounded.  Thus,  we  would  have  been  con¬ 
tent  to  assume,  from  the  outset,  that  H  n  Q  is  bounded. 

We  have  not  done  so,  since  all  of  our  results  (except, 
of  course.  Theorem  9.9)  are  valid  without  this  assumption, 
and,  in  fact,  the  proofs  would  not  be  simpler  were  we  to 
make  this  assumption. 

10.  RELATION  TO  THE  MASS  ACTION  LAWS 

We  now  wish  to  relate  D  to  the  mass  action  laws. 

Recall  from  Secs.  3  and  4  that  the  mass  balance  laws  make 

^For  a  discussion  of  this  problem  for  singlfj-phase 
systems  (chat  is,  for  the  case  K  >=  1)  see  [16j. 
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assertions  about  certain  points  ^  e  E  '  which  satisfy  the 
stoichiometric  conditions.  Since  the  stoichiometric  con 
ditions  are  obtainable  from  the  mass  balance  laws  by 
setting  the  right-hand  sides  of  the  mass  balance  laws 
to  zero,  the  set,  H^,  of  solutions  to  the  stoichio¬ 
metric  conditions  is  a  linear  space;  and  H,  the  mass 
balance  manifold  is  a  translation  of  Hq.  We  accordingly 
define: 

DEFINITION  10.1.  Hq  is  the  set  of  all  points 
of  the  form  M  -  N  where  M  e  H  and  N  €  H.  the 

set  of  all  points  N  e  H,  such  that  N  >  0,  and  such 
that  for  all  9  e  Hq,  we  have: 


I  1  'ki 

k=l  1=1 


+  log 


°k(«y 


=  0  . 


Thus,  Dq  is  the  set  of  points  satisfying  the  mass 
balance,  as  well  as  the  mass  action  laws.  We  have  in¬ 
sisted  that  N  >  0  for  N  Dq,  because  the  mass  action 
laws  as  ordinarily  stated  are  only  applicable  when  N  >  0 
We  now  present  a  theorem  which  generalizes  the  mass 
action  laws.  That  is,  the  condition  of  Theorem  10.2 


reduces  to  the  mass  action  laws  when  N  >  0. 
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THEOREM  10.2.  (Generalization  of  mass  action 
laws.)  N  e  H  n  Q  is  an  element  of  D  if  and  only  if 
Fg(N)  2:  0  for  all  9  for  which  9  is  admissible 

at  N,  where  Fg(N)  is  given  by  Theorem  8.10  or 
Theorem  8.11. 

Proof.  The  theorem  follows  from  the  definition  of 
D  and  the  convexity  of  F.  Note  that  Theorems  8.10  and 
8.11  assure  the  existence  of  Fg(N)  for  0  admissible  at  N. 

THEOREM  10.3.  is  the  set  of  all  N  e  D  with 
N  >  0. 

Proof.  Let  N  e  H,  N  >  0.  We  need  to  show  that  N  e  D 

if  and  only  if  N  e  D^.  Since  N  >  0,  all  9  e  E^  are 

admissible  at  N.  Since  (by  Theorem  8.12)  Fq(N)  *=»  -Fg(N), 

Theorem  10.2  says  that  N  e  D  if  and  only  if  F'(N)  =  0 

y 

for  all  9  e  Hq.  But,  N  >  0  means  that  the  hypotheses 
of  Theorem  8.11  are  satisfied  and  that  the  second  sum 
of  Theorem  8.11  is  zero.  Hence, 


K 


"  I  I 


k-1  j=l 


+  log 


Q.E.D. 


THEOREM  10.4.  If  D  contains  at  .east  one  point 


with  N  >  0,  then  D  is  the  topological  closure  of  D^. 
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Proof.  Apparent  from  Theorem  10.3  and  the  convexity 
cf  D.  Q.E.D. 

11.  THE  POSSIBLE  MINIMUM  SETS 

Thus  far,  I,  K,  and  1^,...,!^^  have  been  fixed  positive 
integers,  H  has  been  a  fixed  linear  manifold  of  E^,  C 
has  been  a  fixed  point  in  E^,  and  hence  D  has  been  a 
fixed  subset  of  E^.  For  the  purpose  of  this  section  and 
this  section  only,  we  will  abandon  the  convention  as  to 
C.  That  is,  we  will  allow  C  to  vary  and  ask  what  possible 
sets  D  may  result.  Our  answer  (Theorem  11.1)  is  that  D 
may  be  any  set  which  is  consistent  with  Theorem  9.2  and 
Theorem  9.8. 

THEOREM  11.1.  Let  N  be  a  fixed  element  of 
H  n  Q.  Suppose  that  N  satisfies  the  condition  of 
Theorem  9.8;  that  is,  suppose  that  for  any  k,j  for  which 
=  0  and  5^  0  we  must  also  have  =  0  for  all 
M  e  H  n  Q.  Then,  it  follows  that  there  exists 
C  e  E^  such  that  the  resulting  minimum  set,  D,  is 
given  by: 


D  =  A(N)  n  H  n  Q  . 
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Proof.  Let  A  be  the  set  of  all  M  €  Q  such  that 


=  0  for  all  k,j  for  which  4  0.  Let  B  be  the 


set  of  all  M  €  A  such  that 


I 

k,j 

Then  B  is  a  bounded  set.  (Note  that  if  N  >  0,  then  B 
is  empty  and  hence  certainly  bounded.) 

Define 


G(M)  =  ^  log  Mkj  -  ^  \(M)  log  ’^(M)  ■ 

k.j  k 

By  Theorem  8.3  (with  C  =  0) ,  G  is  defined  and  continuous 
for  M  €  Q. 

Since  B  £  A  £  Q  and  B  is  a  bounded  closed  set,  G(M) 
is  bounded  on  B.  Hence,  there  exists  a  p  such  that 
G(M)  >  p  for  all  M  e  B. 

Let  M  e  A  and  suppose  that  M  0.  Let 
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G(M  ) 


> 


A  direct  calculation  shows  that  for  any 


t  >  0,  G(Mt)  =  tG(M)  .  Taking  t  =  this  yields 

G(M*)  =  -j—  G(M).  Hence,  G(M)  =  |M|  G(M*)  >  p]M|  . 

If  M  =  0,  then  G(M)  =0.  We  have  thus  proved  that 
G(M)  ^  p|M|  for  any  M  A  and  that  G(M)  =  p|M|  only  when 
M  =  0. 


Let  us  now  define 


N 

=  -  log  — 
k 


kL 

(N) 


if  N,  .  >  0  , 


p  if  N,  .  =  0  . 


Let  b  e  be  admissible  at  N  and  let  ^  i;  Hq.  The 
hypotheses  of  our  theorem  guarantee  that  the  hypotheses 
of  Theorem  8.11  are  satisfied.  Hence,  defining  M  by 
letting  M,  ,  =  .  if  N,  .  =  0  and  M,  .  =  0  if  N,  ,  ^  0,  and 

Kj  kj  kj  kj 

using  the  definition  of  C,  we  get 


F*(N)  =  -  p]Ml  f  G(M)  . 

Since  M  ?  A,  we  have  that  G(M)  ^  Pl^l  >  hence 

F'.(N)  ■■  0,  Furthermore,  Fl(N)  0  unless  M  =  0;  that  is, 

unless  r,  .  =  0  whenever  N,  ,  =  0.  This  shows  that  N  ^  D 
kj  kj 

(by  Theorem  10.2). 
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Let  M'  D.  Then,  M'  ~  N  (by  Lemma  9.5).  Also,  we 
must  have  F  constant  on  the  line  segment  joining  N  and  M' , 
Hence,  =  0.  Thus,  (M’  -  N)|^j  must  be  zero  when¬ 

ever  =0.  This  means  that  must  vanish  whenever 
vanishes.  This  fact  and  the  fact  that  M'  ~  N  imply 
(Theorem  7.5)  that  M'  (N) . 

Thus,  D  <r  A(N).  Since  D  c  H  ^  Q,  we  have 
D  c  A(N)  ^  H  "  Q.  But,  since  N  e  D,  we  must  have 
(Lemma  9.6)  A(N)  n  H  ^  Q  c  D.  Q.E.D. 

12.  EXAMPLES  OF  SPECIAL  CASES  OF  INTEREST 

In  this  section  we  give  some  applications  of  our 

results  to  special  cases.  V/e  give  them  more  for  their 

illustrative  value  than  for  their  importance. 

THEOREM  12.1.  Suppose  that  N  >  0  for  all 

N  €  D;  then  D  contains  at  most  one  point. 

Proof.  We  have  already  shown  that  D  is  the  inter- 

★ 

section  of  some  linear  manifold,  D  ,  with  Q.  If  D 

•k 

contains  two  distinct  points,  then  D  must  contain  the 
line  L  through  these  points.  But  Q  contains  no  complete 
line;  hence,  L  will  contain  both  points  in  Q  and  points 
outside  of  Q,  Therefore,  L  contains  a  boundary  point 
of  Q;  that  is,  a  point  N  with  N  ^  0,  but  not  N  >  0,  a 
contradiction.  Q.E.D. 
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COROLLARY  12.2.  Suppose  that  H  is  such  that 
0  for  all  k  and  all  N  e  H  ~  Q,  and  suppose  that 
H  ''  Q  contains  at  least  one  element,  M,  with  M  >  0. 
Then  D  contains  at  most  one  point,  and,  furthemuDre , 
if  N  €  D,  then  N  >  0. 

Proof.  Theorems  12.1  and  9.8. 

COROLLARY  12.3.  Let  K  =  1  (that  is,  assume  the 
system  consists  of  only  one  phase) .  Suppose  that  H 
does  not  contain  the  origin  (that  is,  suppose  the 
right-hand  sides  of  the  mass  balance  equations  do 
not  all  vanish).  Then,  D  cannot  contain  more  than  one 
point. 

Proof.  By  Theorem  9.2,  if  D  is  non-empty,  it  is  given 
by  D  =  A(N)  ''  H  "  Q.  But,  according  to  the  definition, 

7.4  of  ■ (N) ,  and  because  K  =  1,  the  set  ' (N)  is  the 
set  of  scalar  multiples  of  N.  If  H  were  to  contain  two 
distinct  elements  of  the  line  (N)  ,  H  would  have  to 
contain  / (N)  which  contains  the  origin.  Q.E.D. 

THEOREM  12.4.  Let  N,  N'  <:  D  and  suppose  that 
for  some  k  neither  K,  =  0  nor  N,'  =0.  Then, 


c^(N)  c,^(N')  • 
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That  is,  the  mole  fractions  yielded  by  any  two  points 

in  D  are  the  same  as  long  as  it  is  possible  to  define 
the  mole  fraction. 

Proof.  We  have  N  ~  N' ,  by  Theoreir,  9.2  (or  directly 
from  Lemra  9.5).  Q.E.D. 
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